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When  testing  an  aircraft  model  in  a  wind  tunnel,  the  drag  of  the  model  itself  requires  only  a  low 
proportion  of  the  wind  tunnel  compressor  power  -  say,  the  blockage  ratio,  which  is  a  few  per  cent.  An 
almost  equal  proportion  is  lost  in  the  return  duct.  And  the  balance,  i.e.,  90%  of  the  electric  power 
consumed,  is  spent  to  recover  the  total  pressure  loss  in  wind  tunnel  diffusers.  Therefore,  it  is  crucially 
important  to  reveal  mechanisms  of  the  diffuser  loss  and  predict  it.  Gas  flow  in  a  divergent  diffuser  is 
essentially  multidimensional,  viscous,  often  nonstationary  and  depends  on  flow  nonuniformity  at  the 
diffuser  entry  station;  see  [1-9].  These  difficulties  did  not  allow  researchers  to  identify  major  proc¬ 
esses  in  a  diffuser  with  an  (almost)  optimum  divergence  angle  and  develop  a  rather  clear,  one¬ 
dimensional  loss  theory  which  could  describe  flow  with  general  initial  nonuniformity.  Crocco  in  [2] 
had  cast  doubt  on  even  the  possibility  of  creating  any  one-dimensional  theory  for  analyzing  diffusers; 
nonetheless,  he  appreciated  one-dimensional  approximations  and  had  succeeded  in  this  area. 

Let  us  assume  the  diffuser  internal  flow  to  be  composed  of  two  areas:  the  core  and  the 
stagnant  layer.  The  latter  may  be,  firstly,  the  boundary  layer  over  diffuser  walls,  and  secondly, 
any  stopped  gas  domains  within  the  flow  (such  as  a  model  wake  in  a  wind  tunnel  diffuser).  We 
can  adopt  the  following  assumptions: 

-  within  each  domain  the  gas  parameter  distribution  is  uniform  over  any  cross  section  in  the 
diffuser; 

-  static  pressure  in  the  stagnant  layer  is  equal  to  static  pressure  in  the  core  and  increasing 
downstream; 

-  the  gas  flow  rate  in  the  stagnant  domain  is  insignificant  in  comparison  with  the  core  gas  flow  rate; 

-  the  mean  speed  in  the  stagnant  layer  is  zero;  and 

-  gas  features  no  friction  and  heat  exchange  with  diffuser  walls. 

The  simplified  physical  model  is  depicted  in  Fig.  1;  it  is  very  similar  to  the  flow  pattern  in  a 
divergent  mixing  chamber  of  an  ejector  when  the  ejection  coefficient  is  zero,  so  the  present 
approach  may  be  named  "an  ejection-based  theory". 
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For  the  gas  portion  between  any  two  cross  sections  1  and  2  we  can  write  the  mass  equa¬ 
tion,  energy  equation,  and  momentum  variation  equation.  For  the  inverse  problem  of  aerody¬ 
namics  the  layer  thickness  variation  along  the  longitudinal  axis  can  be  written  as 


^r  =  Fc{x)  =  H  +  ah\ 
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FC(X)=  A-BX 


v  ; 

(where  FI  is  the  relative  thickness  at  the  initial  section,  I+AH  is  the  thickness  at  infinity,  and  X  is 
the  relative  speed).  With  this  in  mind,  the  differential  equation  of  motion  becomes 
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and  may  be  integrated  to  produce 
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Let  us  introduce  a  new  gas  dynamic  function 
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and  use  the  initial  values  to  derive  the  integration  constant  and  interrelate  geometric  and  gas 
dynamic  parameters  of  the  diffuser  as 
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Here,  the  subscript  "H"  is  to  parameters  at  the  initial  section,  v  is  the  diffuser  total  pressure 
recovery  factor,  F  and  Fcore  are  cross-sectional  areas  of  the  diffuser  and  the  flow  core,  and  q(X) 
and  P(X)  are  usual  gas  dynamics  functions. 

If  the  diffuser  has  no  stagnation  domains  (A=A= 0),  then  the  factor  v  equals  unity,  and  dif¬ 
fuser  flow  is  isentropic.  Figure  2  demonstrates  dependence  of  the  diffuser  total  pressure  recov¬ 
ery  factor  on  the  undimensionalized  cross  section  area  F  at  Mach  number  Mh=0.8  and  the 
maximum  possible  increment  of  the  stagnant  layer  thickness,  AH=03.  It  is  seen  that  diffuser 
total  pressure  recovery  factor  gets  lower  if  FI  and/or  the  diffuser  extension  ratio  F  is  increased. 
Thus,  the  analysis  shows  that  total  pressure  loss  takes  place  in  even  a  diffuser  with  no  wall  fric¬ 
tion  and  no  flow  separation  from  walls.  The  loss  occurs  because  the  diffuser  includes  low-speed 
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areas,  and  the  loss  immediately  depends  on  sizes  of  the 
areas.  Any  undertakings  diminishing  the  stagnant  areas 
(for  example,  by  boundary  layer  suction)  always  im¬ 
prove  diffuser  characteristics. 

To  understand  the  mechanisms  decreasing  total  pres¬ 
sure  in  diffusers  and  increasing  their  entropy,  we  may  again 
address  the  diffuser  flow  schematic  in  Fig.  1  ,b.  The  stagnant 
domain  is  subjected  to  the  longitudinal  gradient  of  static 
pressure  and,  thus,  to  the  force  that  acts  against  the  flow. 
For  the  domain  to  be  in  equilibrium  (so  that  the  gas  is  not 
flowing  into  the  diffuser),  the  core  should  apply  to  the  do¬ 
main  a  force  directed  downstream.  Such  a  force  is  provided 
by  friction  between  the  core  and  stagnant  layer.  The  force 
generates  the  kinetic  energy  dissipation  (at  A H=  0): 
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\FcWdP 


The  relative  dissipation  is  proportional  to  the  stagnant  layer  thickness;  and  computation 
shows  that  the  dissipation  is  almost  equal  to  the  relative  energy  of  total  pressure  recovery,  EB  in 
Fig.  3.  Thus,  all  dissipation  in  the  present  theory  is  associated  with  friction  between  the  core  and 
stagnant  layer;  and  there  are  no  additional  mechanisms  that  could  increase  flow  entropy. 

Most  contributions  to  diffuser  flow  theories  note  that  a  considerable  zone  in  a  diffuser  fea¬ 


tures  an  intense  mixing  of  the  flow  by  vortices  which  appear  near  walls.  Also,  many  authors 
notice  that  the  intense  vortex  formation  is  the  primary  factor  in  total  pressure  loss  within  sub¬ 
sonic  diffusers. 


As  shown  in  Fig.  1  ,b,  the  static  pressure  differential  applied  to  the  elementary  volume 
of  the  stagnant  layer  is  equilibrated  by  the  friction  force  which  is  at  a  distance  of  the  stag¬ 
nant  layer  half-thickness,  thus  forming  a  pair  whose  moment  rotates  the  elementary  volume 
around  an  axis  normal  to  the  drawing.  This  moment  is  the  mechanism  that  turbulizes  dif¬ 
fuser  flow  and  causes  the  dissipation  above. 

Let  us  consider  the  simplest  situation  -  a  flat  diffuser  with  a  con¬ 


stant  relative  thickness  of  the  stagnant  layer.  Determine  the  total  turbu- 
lizing  moment  applied  to  the  stagnant  layer.  The  elementary  volume  (in 
Fig.  1  ,b)  is  subjected  to  the  moment 

dM  =dP-Fc^hc=  j FHFHhuHFdP  =  j FHhHPHM2F2dP 

where  FH  and  /?H  are  the  diffuser  section  area  and  height  at  the  initial 
F 

station,  FI  =  —  is  the  relative  area  of  the  stagnant  layer.  The  total 
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turbulizing  moment  is 
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By  introducing  Mj  =  FH,PH  ~^H>we  write  the  relative  turbulizing  moment  as 


MT=^J-  =  H2  FjF2dP=^H2(l-H)- Xf 
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Of  interest  is  a  particular  situation  with  H=  1/7  (that  is,  where  the  stagnant  layer  occupies 
one-seventh  of  the  diffuser  cross  section),  which  produces  the  denominator  exponent  of  2  and 
makes  it  possible  to  evaluate  the  integral  analytically: 
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Line  1  in  Fig.  4  represents  the  turbulizing  moment  dependence  on  the  diffuser  output  speed  at 
Mach  number  A/h=1;  and  line  2,  at  a  low  value  of  M][.  The  greater  the  diffuser  expansion  ratio,  the 
greater  the  total  turbulizing  moment;  so  vortex  formation  is  intense,  and  the  flow  gets  entirely  tur- 
bulized  before  reaching  the  diffuser  output  station. 

It  is  interesting  to  consider  application  of  the  ejection-based  theory  to  the  case  of  incom¬ 
pressible  fluid.  Assume  again  that  the  relative  thickness  of  the  stagnant  layer  increases  linearly, 
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Then  the  pressure  loss  factor  (or  resistance  coefficient)  is 
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.  It  may  be  written  very  simply  in  terms  of  stagnant 


layer  parameters  and  the  input-to-output  speed  ratio.  In  the  case  of  full 
stop  (W=  0)  the  latter  equation  becomes  C  =  H  +  ^ AH  ;  and  if,  in  addi¬ 
tion,  the  stagnant  layer  thickness  is  constant  along  the  diffuser  (A H=  0), 
then  the  pressure  loss  factor  equals  the  relative  thickness  of  the  stagnant 
F 

layer:  C,  =  H  =  — .  This  result  is  the  foundation  for  understanding  the 
F 

performance  features,  computing  the  drag,  and  designing  a  diffuser  with 
a  nonuniform  entry  speed  distribution  (which  is  the  case  downstream  of 
a  wind  tunnel  test  section).  0.01 

Looking  at  Fig.  5,  it  is  clear  that  the  theoretical  total  pressure  re¬ 
covery  factor  v  compares  well  with  experimental  values.  For  instance, 
data  from  [6]  (by  Bedrzhitsky)  obtained  in  an  installation  with  no 
boundary  layer  at  the  diffuser  entry  are  close  to  the  theory  at  H  =  0. 

Experimental  results  by  Lyzhin  and  Pasova  in  [9]  were  obtained  in  a  Figure  4 
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wind  tunnel  with  no  model  in  the  test  section  upstream 
of  the  diffuser;  these  are  close  to  the  theoretical  line  at 
H=  0.1.  When  the  test  section  is  operated  with  free 
circulation  and  an  exit  ledge,  the  diffuser  loss  factor 
gets  increased  and  corresponds  to  calculation  at 
H=  0.15.  When  the  test  section  includes  a  large  model 
of  an  airplane  (whose  wing  area  is  20%  of  the  test 
section  area)  at  20-degree  incidence,  the  loss  factor 
additionally  increases  and  coincides  with  the  theory  at 
H=  0.25.  Indeed,  the  model  in  the  separated  flow 
feature  a  wake  whose  cross-sectional  area  is  about  0.1 
of  the  test  section  area,  and  the  wake  is  added  to  the 
boundary  layer  and  the  ledge  wake. 

Thus,  the  wind  tunnel  diffuser  loss  factor  may  at 
M  =  0  be  described  by  the  equation 
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in  the  case  of  a  long  diffuser  (W=0).  Transonic  flows  can  be  evaluated  by  using 
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where  A=0.3+H,  B= 0.3//,,  and  H  is  the  ratio  of  the  cross- 

sectional  area  of  the  stagnant  layer  to  that  of  the  diffuser  entry. 

The  major  cause  of  pressure  loss  in  both  the  ejector  and  the  diffuser  performance  theory  is 
collision  of  jets.  The  very  idea  to  describe  the  diffuser  loss  in  terms  of  collision  energy  coeffi¬ 
cients  was  proposed  and  actively  developed  by  Abramovich  and  Idel'chik  in  [1,  8].  The  latter 
may  be  regarded  as  predecessors  of  the  present  theory,  although  diffuser  loss  could  be  compared 
with  pressure  loss  in  a  flow  with  a  ledge  in  only  the  case  of  very  large  expansion  angles  of  dif¬ 
fusers.  In  usual  diffusers  of  wind  tunnels  (such  as  those  considered  in  [10])  the  main  stream 
collides  with  its  own  stopped  (boundary)  layer  which  cannot  overcome  the  significant  positive 
pressure  gradient.  And  it  is  of  no  importance  whether  the  stream  thereafter  continues  moving  at 
a  low  speed  in  the  same  direction  (as  in  a  pre-separation  diffuser)  or  slowly  moves  inversely  (as 
in  the  case  of  micro-separation)  -  the  work  of  friction  and  the  work  of  pressure  loss  for  produc¬ 
ing  the  gradient  are  almost  identical. 
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